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Abstract. We present results from a search for gravitational 
microlensing of stars in the Large Magellanic Cloud by low 
Q\ mass objects in the Galactic Halo. The search uses the CCD 
< | light curves of about 82,000 stars with up to 46 measurements 
per night over a period of 10 months. No light curve exhibits a 
form that is consistent with a microlensing event of maximum 
amplification greater than 1.2. This null result makes it un- 
likely that the Halo is dominated by objects in the mass range 
5 10 _8 M o < M < 5 10 _4 M o . 
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The presence of large quantities of "dark matter" in spiral 
galaxies like our own has been inferred from their flat rotation 
curves (Primack et al. 1988). Though a variety of new weakly 
O |interacting elementary particles have been proposed to make 
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up the dark matter, compact objects are also viable candidates 
if they are in a sufficiently dim form. A possible form would be 
objects too light to burn hydrogen (M < 0.07 - O.1M ) (Carr 
1990). 

We report here results from a search for unseen compact 
objects in the Galactic Halo being performed by our collabo- 
ration "EROS" (Experience de Recherche d'Objets Sombres) 
at the European Southern Observatory at La Silla, Chile. Such 
objects can be detected via the gravitational microlensing ef- 
fect (Paczyhski 1986) which would lead to an apparent tempo- 
rary brightening of stars outside our Galaxy as the unseen ob- 
ject passes near the line of sight. The amplification is given by 
A = (u 2 + 2)/[u(u 2 +4) 1 / 2 ] where u is the undeflected "impact 
parameter" of the light ray with respect to the unseen object in 
units of the "Einstein Radius", R E = (AGM d Lx(l - x)/c 2 ) 1/2 . 
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Here, Md is the deflector mass, L is the observer-source dis- 
tance, and Lx is the observer-deflector distance. EROS mon- 
itors stars in the Large Magellanic Cloud (LMC) with L= 
55 kpc yielding typical values Re ~ 2 10 3 i?Q \J Md/ Mq. For 
the "standard" spherical isothermal Halo model (Primack et 
al. 1988), the rate per monitored star for microlensing events 
with amplifications greater than a threshold amplification At 
corresponding to an impact parameter ut has been calcu- 
lated to be 1.66 10" 6 u T (M /M d ) 1/2 yr" 1 (Griest 1991). The 
model assumes a total Galactic "dark" mass of 4.0 10 n M Q 
within 50 kpc of the galactic center yielding a flat rotation 
curve out to the position of the LMC. The time scale for the 
amplification is the time for a Halo object to move through an 
angle corresponding to its Einstein radius and its average is 
t ~ 75 days^/Afd/M©- The result ing achromatic light curve 
has a characteristic shape and, given the small rate for mi- 
crolensing and the preponderance of intrinsically stable stars, 
the event should be the only significant variation on the curve. 
Because of geometry, the events are uniformly distributed in 
impact parameter at maximum amplification, yielding a distri- 
bution of maximum amplifications that falls rather slowly with 
increasing amplification (dN/dA oc 1/A 2 for A > 1). 

We have previously reported results from our Schmidt- 
plate search for long time-scale microlensing (r > 2 days) 
(Aubourg et al. 1993). Two light curves were found that were 
consistent with the microlensing hypothesis. Candidate events 
have been reported by the MACHO collaboration (Alcock et 
al. 1993, Bennett 1994). The OGLE collaboration (Udalski et 
al. 1993) has reported the observation of the microlensing of 
stars in the Galactic Bulge. 

Here, we report on results (described in detail in (Queinnec 
1994) ) from our search for short time-scale (r < 7 days) mi- 
crolensing. It uses a CCD camera to monitor stars in one field 
(1.1° x 0.4°) of the LMC bar. The exposure time was typically 
10 minutes with up to 46 alternating red and blue images taken 
per night. The program is therefore sensitive to microlensing 
events with time scales larger than 30 minutes. About 45000 
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Fig. 1. A sample colour-magnitude diagram. The magnitudes Be and Re correspond to EROS blue and red niters (Arnaud 1994); these 
magnitudes have not yet been related to standard photometric systems. In the present search for microlensing events only the amplification 
(and thus the magnitude change) is relevant. The stars corresponding to the last five events (see text) are shown as filled stars. They are 
all much brighter than the average. 
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that for the blue. Fig. 2a shows the Monte 



Fig. 2. 

The ratio of the s °f the second and first most significant variations for the red light 

Carlo generated events for = 10 — 6 Mq. The accumulation of events in the lower left is due to events with one and only one significant 
variation, as expected for microlensing of stable stars. (Figures obtained for other values of are very similar to fig. 2a.) Fig. 2b shows 
the data. Most events are in the upper right as expected for intrinsically stable stars with no significant variations and for periodic variable 
stars. The position of the cuts are shown in both figures. 



March 31, 1992 while about 82000 stars were monitored be- 
tween August 21, 1992 and March 31, 1993 and again lor the 
same period in 1993 to 1994. This paper presents the first re- 
sults Irom this program using 1800 measurements Irom 1991- 
92 and 4000 measurements Irom 1992-93. Using the calculated 
rate (Griest 1991), the total number ol microlensing events ol 
amplification greater than 1.3 expected to occur during the 
observing time is about O.1(M /M d ) 1/2 , if the Halo is fully 
comprised of objects of mass Md- We can therefore hope to ob- 
serve one or more events if the dominant component of the Halo 
consists of objects in the range ~ 10" 8 M Q < M d < 10" 2 M Q . 
(Below this range the event durations are often less than the 
exposure times.) 

A full description of the CCD camera and telescope can 
be found in (Arnaud et al. 1994). The camera consists of 16 
buttable 579 x 400 pixel Thomson THX 31157 CCDs. Eleven 
CCD's were active in 1991-92 and 15 in 1992-93. The camera 
is mounted on a 40 cm reflector (F/10) refurbished by us and 
the Observatoire de Haute-Provence. 

The first step of the analysis procedure was to construct 
one reference image for each colour and for each season by com- 
bining 50 images taken with good atmospheric conditions. (Be- 
cause of small differences in the filters used and field observed, 
the 1991-92 and 1992-93 data have been treated separately 
but with the same programs.) From the reference images, a 
star finding algorithm then established a star catalogue. After 
elimination of stars near defective pixels, a final star catalogue 
for each season was then formed by matching stars in the two 
colours. The final catalogue then contains about 82000 (45000) 
stars for the 1992-93 (1991-92) season. The limiting magnitude 
corresponds to my ~ 19 with about 30% main sequence and 
70% giants. A sample colour-magnitude diagram from 1 CCD 
is shown in figure 1. 

After the establishment of the star catalogue, light curves 
were constructed by treating individual images. First, each im- 
age was geometrically aligned with the reference using bright 
isolated stars. The star positions in the reference image then 
served as input to a photometric fitting program to determine 
the luminous flux of each catalogue star on the new image. 
Successive images then add one point to the blue or red light 
curve of each star in the catalogue. Photometric errors asso- 
ciated with each point on the curve are estimated empirically 
from the point-to-point variations on a given curve and from 
the overall image quality. They are typically 6 % r.m.s. 

After the elimination of images of poor quality, each 
light curve is subjected to a series of cuts chosen to isolate 
microlensing-like events. As explained below, the efficiency of 
these cuts to accept real microlensing events is determined by 
applying the same cuts to Monte Carlo microlensing events 
that are constructed by amplifying points on randomly selected 
experimental light curves. 

Because of the large volume of data, the first series of cuts 
uses only quantities that can be rapidly calculated. For each 
set of four neighboring measurements in a given colour, we 
calculated a quantity, x|, related to the deviation of the mea- 
sured fluxes from the reference flux, 4> T ef, defined as the most 
probable value on the light curve : 

1=1 



estimated uncertainty. The cuts use the \\ f° r the first and sec- 
ond most significant variations in each colour. Additionally, we 
use the quantity \ 2 calculated as in the above formula except 
that the sum runs over all points not near the most significant 
variation. 

The great majority of stars exhibit only random fluctua- 
tions due to measurement errors. These stars are mostly elim- 
inated by the loose requirement that the most significant vari- 
ation in the blue be within 15 days of the most significant 
variation in the red. Intrinsically variable stars with a very 
significant second variation are eliminated by requiring that 
x\ < 80 for the second variation in each colour. Variable stars 
with long term variations are eliminated by requiring in each 
colour that the \ 2 be less than 2.5 times the number of points. 

After these very loose cuts we are left with about 15 % of 
the original light curves. Because the errors on individual pho- 
tometric measurements are determined only approximately, 
our next cut uses only the ratio of the \\ values for the most 
significant and the second most significant variations. Figure 2 
shows the ratio of the x\ values of the second and first most 
significant variations for the red vs. the same ratio for the blue. 
The Monte Carlo generated events (Fig. 2a) are accumulated 
in the lower left because they have one and only one significant 
variation, as expected for microlensing of stable stars. Fig. 2b 
shows the data. Most events show comparable first and second 
variations as expected for stable stars with variations coming 
only from measurement errors. Requiring that the ratio be less 
than 0.5 in both colours leaves us with only 88 stars. Their 
light curves are then examined in detail and fitted for the the- 
oretical microlensing light curve, neglecting possible star size 
effects. 

Most of the 88 stars show an "unphysical" discontinuous 
flux variation, generally due to inaccurate photometry due to 
bad atmospheric conditions or inaccurate telescope guiding. 
These stars are eliminated by requiring a good agreement be- 
tween the time of maximum variation in the red, Ir, and that 
in the blue, ts- Specifically, we require ts — in < 4<Tt where 
at = 0.05 day \J r / (1 day) is the mean uncertainty in the time 
of maximum of the light curve. After this cut 11 stars remain. 

Six of the remaining stars have variations on long time- 
scales (r > 7days) and are concentrated in regions of the 
colour-magnitude diagram known to contain many variable 
stars. These curves will be discussed in a later publication. 
For the purposes of this paper on short time-scale microlens- 
ing, we make a cut requiring r < 7 days leaving us with five 
stars. This significantly reduces our efficiency for microlensing 
events only if the lensing objects have M d > 1O- 3 M . 

The five remaining stars are indicated in fig. 1; they show 
very small flux variations, of an amplitude comparable with the 
photometric resolution. All events have reconstructed amplifi- 
cations less than 1.16 which, if they were indeed microlensing 
events, would correspond to impact parameters, u > 1.4. Fig- 
ure 3 shows the distribution of fitted impact parameters, u , 
for Monte Carlo events and for the five observed events. In 
contrast to the observed events, the expected distribution for 
microlensing events is concentrated at small impact parame- 
ters. We therefore make a final cut requiring impact parameters 
u < 1.3 leaving no candidates. 

The efficiency of the cuts to accept real microlensing events 
is estimated with Monte Carlo generated lensing events, super- 
imposed on a random sample of the experimental light curves. 



Table 1. The expected number of microlensing events as a 
function of deflector mass, for a standard spherical isothermal 
Halo comprised only of objects of that mass. 
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Fig. 3. The expected distributions of reconstructed impact param- 
eters, u, for M d = 1O _5 M (solid line) and M d = 10~ 7 Mq (dashed 
line). The normalisation of the two distributions is arbitrary, but 
their relative normalisation is correct. The dip observed at small u 
(large amplification) for the distribution at M d = 1O _7 M is due 
to the fact that we fit the theoretical light curve for negligible star 
size, while the Monte Carlo events are generated taking into account 
the actual radius of the source stars. This only results in an over- 
estimation of the fitted impact parameter u but has a small effect 
on the detection efficiency. Also shown are the five observed events 
(hatched area). They are concentrated at higher impact parameters 
(low amplifications). 



The generated events follow the flat geometric distribution of 
impact parameters, u. The relation between u and amplifica- 
tion was modified to take into account two effects. First, the 
finite size of the observed star means that all points on the star 
will not be amplified by the same factor. This effect is impor- 
tant only for M d < 1O -6 M where Re is less than the typ- 
ical stellar radius. The calculated number of expected events 
is reduced by 25 percent for M d = 1O" 7 M . Second, in these 
very dense star fields, stars may be "blended" so that a light 
curve may receive significant contributions from more than one 
star. While this means that we effectively monitor more than 
one star with each light curve, the amplification of a star by 
a given amount A rea i will be reconstructed as a smaller am- 
plification A rec by the photometric programs. This effect has 
been estimated by treating Monte Carlo fabricated images that 
use as input the measured star population in the LMC down 
to luminosities a factor 10 dimmer than the dimmest recon- 
structed by EROS. For our sample of stars, it was found that 
(A rec — 1) ~ P(A rea i — 1) with fi = 0.75 for the brightest star 
associated with the light curve and fi = 0.15 for the second 
brightest star associated with the light curve. The overall ef- 
fect is to reduce the number expected by about 8 (20) percent 
for M d = 1O _6 M (1O- 7 M ). 

Table 1 shows the expected number of events as a function 
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of the deflector mass for a standard spherical isothermal Halo 
comprised only of objects of that mass. The expected number 
of events is greater than 2.3 for 5 10" 8 < M d /M Q < 7 10" 4 
so we exclude this mass range at the 90% C.L. under the as- 
sumption that all objects in the Halo have the same mass. The 
expected number of events is greater than 6.9 for 3 10 - < 
M d /MQ < 1.5 10 - so in this mass range we exclude the pos- 
sibility that such objects could account for as much as one third 
of the Halo. The excluded range applies to any distribution of 
mass that is sufficiently concentrated in the above range. For 
example, we consider a deflector mass distribution of the form 
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and dN/dM = otherwise. Figure 4 shows the excluded 
zone of the parameter space (a, M mtn ). For a > 2 the Halo 
mass is dominated by objects of mass near M mtn and we rule 
out, for a > 3, the range 5 10" 8 < M mtn /M Q < 5 10" 4 . Near 
a = 2 where each decade of mass contains the same total mass, 
the region 10" 12 < M mtn /M Q < 10" 5 is ruled out. For a < 2 
the Halo mass is dominated by high-mass objects and we derive 
no interesting limits. 

The numbers in Table 1 have been obtained for the as- 
sumption of a spherical Halo. The precise mass limits depend 
on the assumed phase space distribution of lenses. Using in- 
stead a flattened Halo (down to c/a = 1/3) does not change 
these numbers by more than 20 percent. Other possibilities 
are discussed e.g. in (Sackett 1993, Giudice 1993, Gould 1994, 
Frieman 1994, Evans 1994). 

In summary, we have searched for microlensing events with 
time scales ranging from 30 minutes to 7 days. The lack of 
candidates in this range places significant constraints on any 
model of the Halo that relies on objects in the range 5 10 - < 
M/Mq < 5 10 -4 . We note that hydrogenous objects of masses 
below this range would have been expected to evaporate before 
the present epoch (de Rujula et al. 1992). We are continuing 
to take and analyze data and expect to improve these results 
soon. 

Acknowledgements. We thank C. Alcock, D. Bennett, A. Bi- 
jaoui, and C. Stubbs for discussions. We thank Ph. Gitton, 
J. -P. Kneib, and the staff of the European Southern Observa- 
tory for help with the observations. 

References 

Alcock C. et al. 1993, Nat 365, 621-623. 

Arnaud M. et al. 1994, Exp Ast 4, 265-278 and 279-296. 



-12 -10 -8 -6 -4 -2 

Log 10 (M mi „ / M sun ) 



Fig. 4. The parameter space (a,-M min ) as explained in the text. 
Inside the outer contour the number of events expected would be 
greater than 2.3 and is excluded at 90% C.L. Inside the inner contour 
the number of events expected would be greater than 6.9 and we 
limit compact objects to less than 1/3 of the total Halo. 

Aubourg E. et al. 1993, Nat 365, 623-625. 

Bennett D. 1994 in Particle Astrophysics, Atomic Physics and 
Gravitation (to be published by Editions Frontieres). 

Carr B. J. 1990, Comments Astrophys. 14, 257-280. 

de Riijula A., Jetzer Ph. and Masso E. 1992, A&A 254, 99-104. 

Evans N.W. and Jijina J. 1994, MNRAS 267, L21-L25. 

Frieman J. and Scoccimarro R. 1994, ApJ 431, L23-L26. 

Giudice G., Mollerach S. and Roulet E. 1993, preprint CERN- 
TH 7127/93. 

Gould A., Miralda-Escude J. and Bahcall J.N. 1994, ApJ 423, 

L105-L108. 
Griest K. 1991, ApJ 366, 412-421. 
Paczyhski B. 1986, ApJ 304, 1-5. 

Primack J.R., Seckel D. and Sadoulet B. 1988, Ann. Rev. Nucl. 

Sci. 38, 751-807. 
Queinnec F. 1994, Ph.D. thesis, Saclay report DAPNIA/SPP 

94-21. 

Sackett P.D. and Gould A. 1993, ApJ 419, 648-657. 
Udalski A. et al. 1993, Act. Astr. 43, 283-294. 



This article was processed by the author using Springer- VerlagI4Tf^X 
A&A style file 1990. 



